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ABSTRACT 

Using B, R, and Ha images of roughly equal-sized samples of low surface brightness 
(LSB) and high surface brightness (HSB) galaxies (~ 40 galaxies apiece), we have 
explored the dependence of Hii region properties on local and global disc surface 
brightness. We have done this by constructing co-added Hii region luminosity func- 
tions (LFs) according to local and central disc surface brightness and fitting Schcchter 
functions to these LFs. The results show that the shape of the Hii region LF within 
LSB galaxies does not change noticeably as different limiting (i.e., fj, > ^um) local 
surface brightness values are used. However, the LFs for HSB galaxies have larger 
values of and are less steep at the faint-end than those of LSB galaxies for limiting 
B-band local surface brightness values as faint as HB,iim —23-24. Both the LFs and the 
data for individual Hii regions show that luminous (L > 10'^^ ergs s^^) Hii regions are 
much more common within HSB discs than within LSB discs, implying that the newly 
formed star clusters are also larger. Taking this into account along with the results 
of Monte Carlo simulations, the shapes of the LFs imply that the regions within LSB 
discs and those within the LSB areas of HSB discs are relatively old 5 Myr) while 
the regions within HSB discs for /xb$24 are significantly younger (< 1 Myr). Since 
the majority of the LSB galaxies do not have noticeable spiral arms and the majority 
of the HSB galaxies do, this may indicate a transition within HSB discs from spiral 
arm-driven star formation to a more locally driven, possibly sporadic form of star 
formation at /is ~ 24, a transition that does not appear to occur within LSB discs. 
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1 INTRODUCTION 

The completeness of the pioneering optical surveys of galax- 
ies suffered from the inevitable bias toward higher surface 
brightness produced by the intensity of the night sky. Dig- 
ital detectors and data reduction techniques used in mod- 
ern times have weakened the effect of this bias. Because of 
this, the study of so called low surface brightness (LSB) 
galaxies has become its own sub-field within extragalactic 
astronomy. Results from this study have since shown LSB 
galaxies to be larger at the sa me luminosity as hig h sur- 
face brightness (HSB) galaxies (|de Blok et al. I [19961 ) while 
still obeying the same TuUy-Fisher relationship, but only if 
the entire baryonic content of the galaxies is taken into ac- 
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count (iMcGaugh et al. Il2000f). LSB discs also tend to be less 
dense (|de Blok fc van de r Hulst 1998) and have higher ga s 
mass fractions than HSB discs (,McGaugh fc de Blok Ill997l ). 
In addition, few LSB galaxies have detect able amounts 
of molecular gas, desp ite several se a rches (lO'Neil et al. I 



2000, 2003; Matthews fc Gao I 1200 ll : lO'Neil et al. I 12004 



[Matthews et al. II2005I ') . LS B galaxies also ten d to have lower 
metallicities (~ 1/3 solar; iMcGaugh Ill994l ) and are more 
frequently dwarf and/or irregular galaxies without notice - 
able spiral arms than HSB galaxies (|Helmboldt et al~ll2004l ). 
In addition to this, the Hi column densities of LSB galax- 
ies tend to be relat ively low, at or below the critical limi t 
for star formation (|de Blok et al. Ill996l : IXennicutt Ill989l '). 
In light of this, it is not surprising that the star forma- 
tion rates (SFRs) for LSB galaxies are typically about ten 
times lower than those of HSB galaxies (jMihos et al. 1 ll999l : 
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iHelmboldt et al. II2004I ). Despite this, LSB discs tend to be 
relatively blue and some have re latively large Ha emiss ion 
line equivalent widths (see, e.g., IHelmboldt et aini2004l ). 

These properties of LSB galaxies imply that despite 
conditions that are not conducive to star formation, many 
LSB galaxies are actively forming stars. Because of this, 
it is possible that they do not form stars in the same 
way as typical spiral galaxies. In fact, it has recently been 
shown that LSB galaxies tend to form fewer high lumi- 
nosity Hll regions and have larger relativ e amounts of dif- 
fuse ionised gas (DIG) than HSB galaxies (IHelmboldt et al~l 
I2OO5I : lOev et al. feoOTi : lO'Neil et al. 2007) . While the larger 
amounts of DIG may reflect differences in the conditions 
of the interstellar medium (ISM) within LSB galaxies, the 
fact that the distribution of Hll region luminosities among 
LSB galaxies is markedly different from that for HSB galax- 
ies points to a possible fundamental difference between how 
star clusters are formed within these two types of galaxies. 
It is remarkable that despite the differences in their star for- 
mation histories, LSB galaxies still manage to form roughly 
the same number of stars at a given circular velocity as HSB 
galaxies. 

From recent results, it is not clear if this difference in 
star formation history results from global or local mecha- 
nisms. In other words, does the lack of spiral arms cause 
star clusters to form differently throughout a typical LSB 
disc than throughout an HSB disc, or do star clusters form 
in a similar way in the LSB regions of HSB discs as they 
do throughout LSB discs? To address this question, we have 
used previously published and newly processed data to com- 
pile two roughly equal-sized samples of HSB and LSB galax- 
ies with B, R, and Ha images. We present here a description 
of the old and new data (§2), the results of a detailed anal- 
ysis of the Hll region properties within these two groups 
of galaxies and how they depend on both local and global 
properties (§3), and Monte Carlo simulations constructed to 
help interpret these results (§4). 



2 GALAXY SELECTION AND 
OBSERVATIONS 

To compare properties of Hll regions within HSB and LSB 
galaxies, we have compiled broad-band B and R images and 
narrow-band Ha images of 89 nearby (14 ^6000 km s~^) disc 
galaxies, 45 HSB galaxies and 44 LSB galaxies. Based on the 
so-called Free man value of /i o.a = 21.65 ± 0.3 for "normal" 
disc galaxies ([Freeman Ill970l ), we have adopted a definition 
for HSB galaxies of fio,B < 21.7. Our working definition of 
an LSB disc galax;y is any galaxy with a B-band central disc 
surface brightness > 22 mag arcsec"'^. This is one magnitude 
brighter than what is typically used (see, e.g., iBothun et al. I 
Il997l : llmpev fc Bothun Ill997h . however, it is still more than 
la fainter than the Freeman value. In addition, as several au- 
thors have demonstrated, the distribution of surface bright- 
ness am ong galax ies is con tinuous a nd not bimodal (e.g ., 
IJansen [[2000 : Blanton et al. 2003: H elmboldt et al. \\2004 ). 
implying this definition is somewhat arbitrary. We there- 
fore chose a limit significantly lower than the Freeman value 
for normal disc galaxies that also yielded a sample of LSB 
galaxies similar in number to our sample of HSB galaxies. 



2.1 Previously published data 

The images were obtained as part of two separate observ- 
ing runs in May and October of 2000 conducted at Kitt 
Peak National Observatory (KPNO) and Cerro Tololo In- 
teramerican Observatory (CTIO), respectively. The CTIO 
run consisted of imaging 69 gala xies selected from t he Hi 
Parkes All Sky Survey (HiPASS: [Barnes et al. II2OOII ) with 
Vr < 2500 km s"\ All 45 HSB galaxies were obtained from 
this sample as well as 20 of the LSB galaxies. The reduction 
and analysis of the CT IO d ata was presented i n det ail by 
[Helmboldt et al. I ([2004[ ) and [Helmboldt et al. I l[2005h , and 
we will briefly discuss them here. The images were calibrated 
using observations of standard stars throughout each night. 
For the narrow-band observations during the last two nights 
of the run, this calibration scheme did not work and a rela- 
tion was derived using the calibration solutions for the other 
Ha and R-band observations given by the following 



ArriHa = ArriR + 2.51og ^ 



Ca— 1 +0.5 

THc 



(1) 



where Am is the difference between the instrumental mag- 
nitude and the actual magnitude, r is the exposure time in a 
particular filter, a is the scale factor each R-band image was 
multiplied by to be used to subtract the continuum emission 
from the Ha image, and C is a constant empirically deter- 
mined to be 0.73 ± 0.15. Here, we define the instrumental 
magnitude to be — 2.51og DN/r -I- 25 where DN is the flux 
measured from the image in ADUs. 

The 25 mag arcsec"'^ isophote was then identified on 
each calibrated B-band image and an ellipse was fit to it, ex- 
cept for one galaxy for which the 26 mag arcsec"^ isophote 
was used because it had /xo.s ~ 25. The position, ellipticity, 
and position angle of the fitted ellipse was then used to con- 
struct elliptical apertures within which the median B and 
R band surface brightnesses were measured to estimate the 
surface brightness profile of each galaxy. The B and R pro- 
files were then simultaneously fit with an exponential disc 
model for all radii where the surface brightness was between 
0.3 and 2.472 mag above the la limiting isophote. This mag- 
nitude range was chosen because for an exponential disc, this 
corresponds to a span of two disc scale lengths. Each fit was 
performed with a standard weighted linear least squares rou- 
tine to obtain the central surface brightness values in each 
band and a single value for the disc scale length, h, as well 
as error estimates for all three quantities. The fits were then 
used to obtain total extrapolated fluxes by integrating each 
exponential fit from the radius at which the surface bright- 
ness is 2.172 mag above the la limiting isophote to infinity, 
and then adding this fiux to the total flux within that radius 
measured from the image. 

Continuum subtracted Ha images were made by scal- 
ing the R-band images so that the fluxes of several stars 
on the R-band and Ha images matched, on average, and 
subtracting the scaled R-band image from the Ha image. 
The Ha fllter transmission curve and the galaxy redshift 
were then used to determine the Ha-|-[NII] flux from each 
calibrated, continuum subtr acted Ha im age. Using the emis- 
sion line measurements o fTremonti et al. (200 4) from Sloan 
Digital Sky Survey fSDSS: [York et al. [I2OO0I ) galaxy spec- 



tra, we have derived relationships between B— R colour in- 
dex and the ratios of the Ha flux and equivalent width 
to the observed Ha+[NII] flux and equivalent width. The 
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iTremonti et al. I (|2004l ') measurements included a linear com- 
bination of model spectra fit to each galaxy's spectrum so 
that the effects of stellar absorption would be fully accounted 
for. We also corrected both th e Hq fl ux and equivalent width 
measured bv lTrernonti et ah" ||2004|) for internal dust extinc- 
tion according to ICalzetti I ( 200lh using the ratio of Ha to 
H/3 emission line flux to compute E(B— V) for the ionised 
gas, assummg_ Ha/H/3 = 2.87 (for T= lO'^ K and case B re- 



combination; 



Qsterbrock &: Ferland I [2OO6I ) . Using emission 
line ratios, we identified 12,262 ne arby (m^ < 16) activel; 



star forming galaxies acc ording to Helmboldt et al. I (|2008i ) 
from the SDSS data of ITremonti et al. I (|2004l l and have 
measured B— R indices from their spectra. In Fig. [T] we 
have plotted the ratios of the Ha flux and equivalent width 
to the observed Ha-|-[NII] flux and equivalent width versus 
B— R for these galaxies with median values plotted in bins of 
B— R. Also plotted are broken power-law functions fitted to 
the data (not the median values) which have the following 
form 

X = (0.42 ±0.066)(B-i?) + (-0.28 ±0.085) 
if {B-R)> 0.992, 
= (0.17 ±0.11)(B-i?)-f (-0.032 ±0.11) 

if 0.992 > (B - i?) > 0.25, 
= (0.011 ±0.11) if (B-i?)< 0.25. (2) 
Y = (0.16 ±0.066)(B-i?)± (-0.18 ±0.084) 
if (B - ii) > 1, 
= (-0.020 ± 0.086) if (B - 7?) < 1. (3) 

where X = log F{Ha)o/F{Ha + [NII])obs and Y = 
log W{Ha)o/W{Ha + [NII])obs- We must note that the 
SDSS spectra were obtained using a flbre-fed spectrograph, 
with each flbre having a 3 arcsec aperture. This implies that 
many of the spectra may be bulge-dominated and the above 
corrections may be more appropriately applied to nuclear re- 
gions of galaxies. To check this, we obtained Se r sic in dices, 
Us, for the SDSS galaxies from *Bla nton et al. I (|2003l ) and 
re-computed the median values plotted in Fig. [1] for galax- 
ies with 0.5 < ris < 1.5, i.e., likely disc-dominated galaxies. 
These new median values are plotted in Fig. [1] as large, open 
red circles and are not significantly different from the me- 
dian values computed for the entire sample of star forming 
galaxies. We therefore conclude that it is reasonable to ap- 
ply the corrections given in equations (2) and (3) to galaxy 
discs based on their B— R colours and not on integrated or 
nuclear colour indices. 

We found no significant correlation between the ratios 
plotted in Fig. [1] and either absolute magnitude or surface 
brightness. We also note that different e mpirically derived 
correc tions were used for th e se da ta by iHelmboldt et alH 
l|2004h and [Helm bol dt et al. I (|2005l ). However, those previ- 
ous relations did not include the effects of stellar absorp- 
tion, and we have therefore concluded that equations (2) 
and (3) provide better, more extensive corrections for the 
Ha±[NH] data. These new empirical relations with B— R 
colour were used to convert the measured Ha±[Nn] fluxes 
to extinction corrected Ha fluxes using the B— R colour of 
each galaxy's disc computed using the exponential fits to 
the B and R surface brightness profiles described above (i.e., 
(B - R)disc = Mo,s - /io.ifO- We used the {B - R)disc colours 
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Figure 1. The ratio of the Hq flux (upper) and equivalent width 
(lower) to the observed Ha-|-[N II] flux for 12,2 6 2 star forming 
SDSS galaxies as measured by ITremonti et al. includ- 
ing co rrections for internal dust extinction according to lCalzetti I 
1I2OOII) . versus B— R colour index as measured from the SDSS 
spectra (see §2.1). Also plotted are median values within bins of 
B— R and broken power-law fits to the data [not to the median 
values; see equation(2)]. Median values for galaxies with Sersic 
indices from lBlanton et al. ] Hooi) between 0.5 and 1.5 (i.e., disc- 
dominated galaxies) are also plotted as large, open red circles. 



because (1) the vast majority of the star formation within 
these galaxies occurs within their discs, and (2) the expo- 
nential fits yield B— R colour estimates for the discs that 
are reasonably free from the effects of any substantial small 
scale colour fluctuations and bulge or bar contamination. 

To measure the properties of the Hll r e gions , the 
Hllphot algorithm developed by iThilker et al. I l|2000l) was 
used. The algorithm is designed to identify Hll regions by 
smoothing the Ha images with kernels of different sizes and 
fltting models to signiflcant peaks in the smoothed images. 
The boundary of each Hll region is then grown until a limit- 
ing emission measure gradient is reached (a limit of 1.5 EM 
pc~^ was used), or until the boundary reaches the boundary 
of another region. Following this, surface flts to designated 
background pixels are used to estimate the amount of dif- 
fuse emission a long the lines of sight passing through each 
Hll region (see iThilker et a"ni2000l ; IHelmboldt et al. II2005I . 
for more detailed descriptions). 



2.2 Newly processed data 

To supplement the imaging data for LSB galaxies avail- 
able from the CTIO observing run, a sample of 32 
galaxies previously identified as LSB disc galaxies was 
observed at KPNQ. The 32 targets we r e culled from 
the l it erature (ISchom mer fc Both un ' 'l983'; 'Bothun et ah_ 
19851 ; iKraan-K ortcwcg l Il986l; Elmcgrccn & Elmcgrccn I 
19871 ; ITuUv fc Fisher I Il988l ; ISchombert fc Bothun I ligSsT 
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Figure 2. In the upper panel, the total extrapolated B-band 
magnitude (see §2.1) for the 12 LSB galaxies observed at KPNO 
that are in the RC3 catalogue versus the RC3 B-band magnitudes. 
The difference between the B-band magnitude measured using 
our KPNO images and the RC3 magnitude is plotted versus the 
RC3 magnitude in the lower panel. In each panel, the dashed line 
represents the case where the magnitudes are the same and is not 
a fit to the data. 

ISchneider et al. 1 ll990l : IO'Neill et al. |[l997l ) and were chosen 
to represent the wide variety of known LSB disc galaxies, 
ranging from dwarf (disc scale lengths ~500 pc) to giant 
(scale lengths ~10 kpc) LSB galaxies. However, because of 
cloudy conditions during the observing run, we did not per- 
form standard star observations and the usual method of 
calibrating the images could not be used. In the interven- 
ing years, the completion of the bulk of the imaging for 
the SDSS has provided a means to calibrate the B and R 
band images of 24 of the LSB galaxies observed. We have 
done this by measuring the fluxes of several (~6-10) stars 
on the KPNO B and R images which have measured g and 
r magnitudes from SDSS images. The g and r magnitudes 
were con verted to B and R m agnitudes using the conversions 
given by ISmith etHI (|2002D . and the mean difference be- 
tween the raw KPNO fluxes and these magnitudes was used 
to calibrate each image. The typical uncertainty in these cal- 
ibration solutions is about 0.06 mag and 0.02 mag for the B 
and R bands, respectively. 

There are 12 galaxies from our KPNO sample with B- 
band magnitudes from the RC3 catalogue. As Fig.[2]demon- 
strates, our total B-band magnitudes (see below) agree with 
the RC3 values to within about la on average and that the 
difference between the two is not systematic and does not 
depend on the brightness of the galaxy. Note that the la 
uncertainties in our total B-band magnitudes are combina- 
tions of the calibration errors, the photon errors, and the 
errors in the parameters from the fits to the surface bright- 
ness profiles used to compute fiuxes integrated to radii of oo 
(see below). 
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Figure 3. The B-band central surface brightness (left) and B— R 
disc colours (right; see §2.1) distributions for the HSB (grey) and 
LSB (light grey) galaxies. The distributions for the LSB galaxies 
from the CTIO data only (black histograms) are also included. 

For consistency, the surface photometry for these new 
images was performed in the same way as the CTIO images, 
with one important difference. Since many of these galaxies 
have relatively low central surface brightnesses, the B-band 
25 mag arcsec"^ isophote was not ideal to use to charac- 
terise the shape of each galaxy if it could be identified at 
all. Instead, a different isophote was used for each galaxy, 
determined by eye through trial and error by varying the 
isophote level in 0.5 mag steps. The isophotes used are listed 
in Table [2] and range from 26 to 27 mag arcsec"'^. We have 
also listed in Table [2] the ellipticity and position angle of 
this isophote along with the exponential disc scale length, 
the B and R central surface brightnesses, and the total ex- 
trapolated B and R magnitudes, all deter mined in the same 
way as the CTIO data (see above and iHelmboldt et alTI 
[2004') . The surface brightness profiles for all 24 galaxies are 
displayed in Fig. |4] along with the exponential disc fits. The 
/io,B distributions plotted in the left panel of Fig. [3] show 
that the inclusion of the newly processed KPNO data not 
only doubles the number of LSB galaxies, it also provides a 
much better sampling of the 23.5 < mo.b < 24.5 region. In 
the right panel of Fig. |3l we have plotted the distributions of 
{B — R)disc- These distributions show that the inclusion of 
the new KPNO data has also introduced a significant num- 
ber of relatively red LSB discs that supplement the relatively 
blue LSB discs found among the CTIO data. These redder 
discs have made the {B — R)disc distributions for the HSB 
and LSB galaxies essentially the same. 

Without standard star observations or spectra, we have 
no direct means of calibrating the Ha images obtained for 
these galaxies. However, for all but five of the galaxies ob- 
served a KPNO, the same set of 30 A wide filters, kindly 
provided by Richard Rand, was employed as was used for 
most of the CTIO observations. For the remaining galaxies. 
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Table 1. Properties of LSB Galaxies Observed at KPNO 



R.A. {J2000) Doc. {J2000) Vr PA h log F(Ha) log W{Ha) 

Name (h m s) (o ' ") (km s ^) e {°) f^iso C) Mo B Mo R (ergs/s/cm'^ ) (A) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



NGC 3447 


10 


53 


26 


,8 


+ 16 


46 


44, 





1067 





,32 


— 73.0 


26. 





27.2 


22 


,5 


21 


,6 


13 


,1 


12. 


3 


— 11, 


9 


1 


,71 


UGC 6151 


11 


05 


56 


,3 


+ 19 


49 


31. 





1335 





,18 


88.0 


27. 





18.2 


23 


,2 


22 


,3 


14 


,8 


13. 


9 


— 13, 





1 


,33 


Ul-4 


11 


38 


25 


,8 


+ 17 


05 


03, 


1 


3452 





,38 


69.9 


26. 





3.4 


22 


,1 


21 


,1 


17 


,3 


16. 


3 


— 13, 


7 


1 


,11 


NGC 4688 


12 


47 


46 


5 


+04 


20 


09. 


9 


986 





,16 


39.7 


26. 


5 


26.6 


22 


,3 


21 


,4 


12 


,9 


12. 


1 


— 11, 


9 


1 


,49 


LSBC F574-10 


12 


50 


31 


,4 


+ 17 


28 


18. 


7 


844 


0, 


,72 


-74.5 


26. 


,0 


6.1 


22 


,9 


22 


,1 


16, 


,8 


16. 


1 


-14, 


,0 


0, 


,94 


UGC 8011 


12 


52 


21 


1 


+21 


37 


46. 


2 


776 


0, 


,08 


-12.2 


27. 


,0 


19.7 


23, 


,9 


23, 


,0 


15, 


,3 


14. 


,4 


-13, 


2 


1 


,43 


UGC 8061 


12 


56 


44, 


,0 


+ 11 


55 


55. 


,3 


562 


0, 


,17 


-86.0 


27. 


,0 


16.6 


23, 


,8 


22 


,5 


15, 


,5 


14. 


,3 


-13, 


5 


1 


,22 


UGC 8091 


12 


58 


40, 


4 


+ 14 


13 


02. 


,9 


214 


0, 


,39 


-74.1 


26. 


5 


16.2 


23, 


,0 


22 


,2 


14, 


,7 


13. 


,9 


-12 


,0 


2 


,16 


UGC 8155 


13 


03 


14, 


7 


+07 


48 


07. 


,9 


2925 


0, 


,24 


-86.2 


27. 


,0 


30.8 


23, 


,5 


22 


,2 


13, 


,6 


12. 


,4 


-12 


5 


1 


,07 


UGC 8441 


13 


25 


29, 


,1 


+ 57 


49 


20. 


,0 


1519 


0, 


,17 


1.6 


27. 


,0 


27.2 


23, 


,6 


22 


,6 


14, 


,3 


13. 


,4 


-12, 


,8 


1 


,31 


NGC 5409 


14 


01 


46, 


,1 


+09 


29 


25. 


,0 


6259 


0, 


,01 


12.1 


26. 


,0 


23.8 


22 


,9 


21 


,5 


13, 


,7 


12. 


,3 










UGC 9024 


14 


06 


40, 


,6 


+22 


04 


12. 


1 


2323 


0, 


,16 


-26.8 


27. 


,0 


13.5 


23, 


,6 


22 


,6 


15, 


,6 


14. 


,6 


-13, 


,3 


1 


.26 


UGC 9245 


14 


25 


26, 


,9 


+ 56 


19 


12. 


,6 


1885 


0, 


,30 


-79.3 


27. 


,0 


23.0 


23, 


,3 


22 


,2 


14, 


,3 


13. 


,3 


-12, 


7 


1 


.22 


UGC 9380 


14 


34 


39, 


,2 


+04 


15 


46. 


,0 


1693 


0, 


,32 


-66.5 


27. 


,0 


20.2 


23, 


,3 


22 


,5 


14, 


,7 


13. 


,9 


-12 


,8 


1 


.50 


NGC 5866B 


15 


12 


07, 


,2 


+55 


47 


06. 


,3 


841 


0, 


,45 


-82.5 


26. 


5 


24.3 


23, 


,0 


21 


,9 


14, 


,0 


12. 


,8 


-12, 


7 


1 


.26 


LSBC F583-2 


15 


37 


00, 


,5 


+ 20 


47 


42. 


,0 


1719 


0, 


,36 


38.8 


27. 


,0 


11.3 


23, 


,9 


22 


.9 


16, 


,4 


15. 


5 


-13, 


7 


1 


.50 


NGC 5964 


15 


37 


36, 


,3 


+05 


58 


24. 


,3 


1447 


0, 


,09 


-27.4 


26. 


,0 


43.2 


22 


,5 


21 


.4 


12 


,0 


11. 


,0 


-12, 


1 


1 


.01 


LSBC F583-5 
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Table 2. Col. (4): Radial velocity. Col. (5): Isophotal ellipticity. Col. (6): Isophotal position angle (from north through east). Col. (7): 
B-band surface brightness of isophote used for Col. (5) and (6). Col. (8): Exponential disc scale-length. Col. (9): Disc central B-band 
surface brightness. Col. (10): Disc central R-band surface brightness. Col. (11): Total extrapolated B-band magnitude (see §2), corrected 
for foreground dust extinction. Col. (12): Total extrapolated R-band magnitude (see §2), corrected for foreground dust extinction. Col. 
(13): Log of the total Ha flux (see §2). Col. (14): Log of the total Ho emission line equivalent width (see §2). 



two broader Ha filters provided at KPNO were used which 
are similar to those available at CTIO which were employed 
for 21 of the 69 galaxies observed there. Because the same or 
similar Ha filters were used and since we were again using 
scaled versions of the R-band images to subtract the con- 
tinuum emission from the Ha images, it was possible to use 
the derived R-band calibrations with equation (1) to cali- 
brate the Hq images. We have done this assuming the mean 
value of C = 0.73 determined using the CTIO observations 
and have also included the corrections for internal dust ex- 
tinction, [NH] contaminat i on, an d stellar absorption derived 
from the lTremonti et al. I l|2004l ) data given in equations (2) 
and (3), again using the disc B— R disc colours computed 
using the exponential fits to the surface brightness profiles 
(see §2.1). We note that the uncertainty in the value of C 
corresponds to an uncertainty of about 0.09 dex, which is 
comparable or less than the uncertainties in the corrections 
given by equations (2) and (3) . The resulting calibrated con- 
tinuum subtracted Ha images were then used to measure the 
total Ha flux and equivalent width for each galaxy, which 
is listed in Table [S] We have also used Hllphot to identify 
and perform photometry on Hll regions for each of these 
images using the same settings as were used for the CTIO 
images. We note that for one of the LSB galaxies observed 
at KPNO, NGC 5409, no Ha image was obtained. The Hll 
region luminosity functions (LPs) of the 23 galaxies with Ha 
images were computed using the Hllphot results, assuming 
Ho =70 km Mpc~^, and are displayed in Pig. [5] LPs for 
galaxies with radial velocities < 1000 km s~^ for which the 
distances may be relatively uncertain are highlighted in red. 



The Hll region properties and LPs of the galaxies observed 
at CTIO are presented in lHelmboldt et al. I (|2005l ). 



3 HII REGIONS 

3.1 Global properties 

With B, R, and Ha images for 45 HSB and 43 LSB galaxies, 
we are in a good position to explore the differences in Hll 
region properties between the two classes of galaxies. In Fig. 
[51 we have plotted the total Ha emission line equivalent 
width as a function of B— R index for the HSB and LSB 
galaxies as well as similar data from the Nearby Pield Galaxy 
burvey (NPGS; [j ansen 1 [2000 ). which includes a wide range 
of morphological types. Prom this plot, it can be seen that 
the amounts of star formation relative to the stellar continua 
within both types of galaxies are similar. Both are bluer on 
average than the NPGS galaxies. Por the most part, both the 
LSB and HSB galaxies have similar Ha equivalent widths at 
the same colour indices as the NPGS galaxies. 

However, we start to see distinct differences between 
HSB and LSB galaxies when we look at Hll region prop- 
erties. In Fig. [T] we have plotted the distributions for two 
properties (1) the Hll region covering factor, the ratio of the 
total area occupied by Hll regions taken from the Hllphot 
results to the total disc area taken to be 'dirh? , where h is 
the disc scale length, and (2) the mean luminosity of the 
three brightest Hll regions, L^. We note that for the su rface 
brightness profiles presented in lHelmboldt et al. I (|2004l ) and 
in Fig. |4j the radii are the geometric mean radii of the ellip- 
tical apertures used, i.e., R = \/ab. Therefore, the assumed 
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Figure 4. The R (open red points) and B (solid blue points) band surface brightness profiles for the LSB galaxies observed at KPNO; 
no offsets have been applied to the data. The best fitting exponential profiles (see §2.1) are plotted in red and blue for the R and B band 
profiles, respectively. The Icr limiting isophote is plotted as a horizontal dashed line for each band. 



total disc area of 9Tvh^ is the area of an ellipse with a semi- 
major axis of 3h{b/a)~^^^ , where b/a is the axis ratio of 
the elliptical apertures used. No attempt has been made to 
de-project the estimated total disc area or the area occu- 
pied by Hll regions. However, we have found no significant 
correlation between b/a and either the Hll region covering 
factor or the average number of regions per square kilopar- 
sec with Spearman rank correlation coefficients of -0.08 and 
-0.18, respectively. Additionally, the mean axis ratios for the 



LSB and HSB samples are similar {<b/a>— 0.65 and 0.60, 
respectively). It is therefore unlikely that this has had a sig- 
nificant effect on the results presented here. 

As Fig. [7] demonstrates, the Hll region covering factor 
is significantly lower on average for LSB galaxies than it is 
for HSB galaxies, and the same is true for L3. These results 
are generally consistent with what has been found previously 
llGreenawalt 11 19981: iHelmboldt et al. II2005I : lOev et al. II2OOT : 
lO'Neil et aUbOO?! ). 
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Figure 5. The Hll region luminosity functions (LFs) for the LSB galaxies observed at KPNO. LPs for galaxies with radial velocities 
<1,000 km s^^ for which the distance estimates are somewhat unreliable are highlighted in red. 



To compare the typical distributions of Hll region lu- 
minosities, we have constructed co-added Hll region LFs for 
both the LSB and HSB galaxies. First, we must note that we 
only have radial velocity distances for our galaxies which im- 
plies that those with smaller radial velocities will have larger 
distance uncertainties. These large errors may artificially al- 
ter the shape of any co-added LF, and we therefore must 
exclud e such galaxies from this process. iHelmboldt et alH 
l|2005ll determined that the optimum radial velocity cut-off 
for a co-added LF with bins of 0.2 dex in width is Vr < 840 



km s~^. For all co-added LFs discussed in this paper (see 
§3.2 as well), we have excluded galaxies with Vr < 840 km 
an d have used the same bin width as IHelmboldt et al. I 
(2005!). There are 40 of our HSB galaxies and 37 of our LSB 
galaxies that have Vr > 840 km s~^ . 

After applying the radial velocity cut, the next step to 
co-adding the galaxies' LFs was to determine a limiting Hll 
region luminosity for each galaxy by fitting a power-law to 
the Hll region luminosity as a function of the signal-to-noise 
ratio, S/N, of the integrated flux of the region as determined 
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Figure 6. The total Ha emission line equivalent width versus the 
B— R colour index for HSB (grey) and LSB (light grey) galaxies; 
LSB galaxies observed at KPNO are plotted as open points. For 
comparison, the galaxies from the Nearby Field Galaxy Survey of 
Jansen (2000) are plotted with X symbols. 

by Hllp/ioi. The limiting luminosity was then taken to be 
the luminosity given by the power-law fit for S/N= 5, plus 
the rms deviation of the data about the power-law fit. We 
then constructed an Hll region LF for each galaxy above this 
limiting luminosity. Each LF was then divided by the area of 
the disc, Qtt/i^, to make sure that larger galaxies with more 
Hll regions would not dominate the shape of the co-added 
LF. As noted above, this estimate of the total disc area has 
not been de-projected to the face-on disc area. However, as 
also noted, there is no significant correlation between the 
isophotal axis ratio, b/a, and the average number of Hll 
regions per square kiloparsec, and the mean values oi b/a 
for the LSB and HSB samples are similar. 

For each galaxy group, the average number of regions 
per square kiloparsec was computed for each luminosity bin 
by averaging over all galaxies with limiting Hll region lumi- 
nosities larger than the upper boundary of the bin. Uncer- 
tainties for each LF were computed using jackknife resam- 
pling. Specifically, the galaxies were put in groups of four 
and the LF was computed 11 different times, each time ex- 
cluding one group of four, and the standard deviation for 
each luminosity bin over these 11 computations was taken 
to be the uncertainty. For instances where the shot noise 
{\/Ntot, where Ntot is the total number of regions among all 
co-added galaxies) was larger than the jackknife value, the 
shot noise was adopted as the uncertainty. 

The HSB and LSB galaxy co-added LFs are plotted 
with these uncertainties in Fig. [S] We can see that while 
the LSB LF is well approximated by a single power-law over 
more than two decades in Ha luminosity, the HSB LF shows 
a signific ant amount of curvature. This is similar to what was 
found bv lHelmboldt et ai~l |200l). To quantify this, we have 
fit to the LFs the following function 



Figure 7. Distributions for the Hll region covering factor (see 
§3.1; upper) and the mean Ho luminosity of the three bright- 
est Hll regions (lower) for the HSB (grey) and LSB (light grey) 
galaxies. 

commonly referred to as a Schechter function where a is the 
faint-end slope and L* characterises the amount of curva- 
ture at the bright-end. The Schechter function was chosen 
because it is basically a truncated power-law which is qual- 
itatively what one would expect for an aged population of 
star clusters with a power-law distribution of initial masses 
(see the Monte Carlo simulations described in §4). The fit- 
ted functions are also plotted in Fig. |8] The values of a and 
L. from these fits are -2.08 and 5.79 x 10^® ergs s"^ for the 
LSB LF and -1.52 and 1.86 x 10^^ ergs s"^ for the HSB LF. 
We note that this difference in LF shape is unlikely to be 
artificially caused by the application of the corrections given 
in equation (2) since these corrections are solely functions of 
{B — R)disc, and Fig.|3]demonstrates that the {B — R)disc dis- 
tributions for the LSB and HSB galaxies are essentially the 
same. We must also acknowledge that the shape of the LF, 
particularly for HSB galaxies for which the Hll region cover- 
ing factor is large r, may be artificially fiattened by blending 
of Hll regions (see lScoviUe et al. II2001I ). We will address this 
issue in more detail in §3.3. 

3.2 Local properties 

While the results presented in Fig. [7] and |S] mostly reinforce 
earlier results, a question still remains regarding whether 
they arise from differences in the global properties of HSB 
and LSB galaxies, or from local conditions within the LSB 
and HSB regions of all discs . While results previously pre- 
sented bv iGreenawalt I (Il998l ) suggest that these properties 
vary with local surface brightness within HSB discs, we now 
have the ability to explore whether or not that variation is 
seen in LSB discs and if it is indeed similar to what has 
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Figure 8. Co-added Hll region luminosity functions (see §3.1) for 
tlie HSB (grey) and LSB (light grey) galaxies with radial velocities 
>840 km s~^ (see §3.1). Schechter functions fit to the LFs are also 
plotted as dashed lines. 

been seen for HSB discs. We have done this by measuring 
the local B and R band surface brightness in the area near 
each Hll region using a circular aperture with a radius of 
2y/n^/-K where Up is the number of pixels occupied by the 
Hll region on the image. We have defined the sizes of the 
apertures in this way rather than using an aperture with a 
fixed physical size to ensure that each aperture was small 
enough to represent the local conditions around each region 
while being large enough that the region itself did not dom- 
inate the local fiux. In Fig. [9l we have plotted the Hll region 
luminosity versus local surface brightness. For both LSB and 
HSB galaxies, there appears to be a trend for Hll regions to 
be more luminous in higher surface brightness areas. The 
fact that there are few low luminosity Hll regions in higher 
surface brightness areas may be a selection effect as they 
are harder to detect in such areas. However, the paucity of 
highly luminous Hll regions in more LSB areas cannot be 
the result of such a bias, implying that while this possible 
selection effect may infiuence the observed trend, it cannot 
be the cause of it. In Fig. |5] for /is < 23, the trend for 
LSB galaxies becomes somewhat steeper than that for HSB 
galaxies. However, this is heavily influenced by a few LSB 
galaxies which contain some HSB regions and this differ- 
ence is not seen in the R-band data. In both the B and R 
data, the trend appears to become steeper for fiB < 22.5 
and fiR < 21.5. 

To further explore and constrain the results shown in 
Fig. |9l we have fit a line to log L{Ha) as a function of 
local surface brightness for each galaxy to obtain a gradi- 
ent. The distributions for these gradients using both the B 
and R data are plotted in Fig. [10] separately for LSB and 
HSB galaxies. On average, the log L{Ha) gradients for LSB 
galaxies appear to be about 0.1 dex mag~^ arcsec'^ less steep 
than those for HSB galaxies with both distributions peaking 



Figure 9. The Ha luminosity of Hll regions from LSB (loft) and 
HSB (right) galaxies as a function of local surface brightness in 
the B (upper) and R (lower) bands. In the left panels, data for Hll 
regions within the LSB galaxies observed at KPNO are plotted 
as open circles. Mean values are plotted within surface brightness 
bins for both LSB (light grey) and HSB (grey) galaxies. 

at values < 0. This is likely a refiection of the fact that the 
LSB galaxies have a substantially higher fraction of regions 
with fiB > 22.5 and fiR > 21.5, where the trends seen in Fig. 
[9] are less steep, than HSB galaxies. We note that the gradi- 
ent distributions for both the B and R bands peak around 
—0.15 to —0.05 dex mag~^ arcsec^, which is similar to the 
slopes for the median values plotted in Fig. [O] 

The results illustrated by Fig. [9] and [10] imply that for 
a given local surface brightness, the typical Hll region lu- 
minosity for LSB galaxies is likely similar to that for HSB 
galaxies. This implies that local surface brightness/density 
strongly affects star cluster for mation within galaxy discs, 
as suggested by the toy model of lO'Neil et al. I (|l998l ) . How- 
ever, this does not mean that global disc properties do not 
have a significant effect on the detailed properties and for- 
mation histories of star clusters. For instance, at the same 
local surface brightness, the shapes of the distributions of 
Hll region luminosities are not necessarily the same for LSB 
and HSB galaxies, and any difference in shape may point to 
the inffuence of global properties on cluster formation. To 
explore this further, we have constructed co-added LFs for 
Hll regions based on local B-band surface brightness sepa- 
rately for HSB and LSB galaxies. We did this by making 
three co-added LFs separately for LSB and HSB galaxies 
using three different limiting local B-band surface bright- 
ness values, fiB.iim =22, 23, and 24 mag arcsec"^. In other 
words, each co-added LF was remade using only Hll regions 
with local fiB > l^B,iim- The six resulting LFs are plotted 
in Fig. [11] along with Schechter function fits. From these 
plots, it can be seen that the HSB LFs have a noticeably 
larger amount of curvature and a faint-end with a shallower 
slope than the LSB LFs for ^B.Um < 24. This is better il- 
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Figure 10. The distributions of the gradients of the log of the 
Hll region Ha luminosity (see §3.2) as a function of local surface 
brightness in the B (upper) and R (lower) bands for Hll regions 
within LSB (light grey) and HSB (grey) galaxies. 



Figure 11. Co-added LPs for Hll regions from HSB galaxies (grey 
to red) and from LSB galaxies (light grey to blue) using different 
limiting local B-band surface brightness values (see panel). 



lustrated in Fig. [TS] where we have plotted the parameters 
of the Schechter functions fitted to the LFs in both Fig. [S] 
and llll as functions of ^B.Um- While the faint-end slope, a, 
remains roughly constant with fiB,iim for LSB galaxies, it 
appears to become steeper with decreasing limiting surface 
brightness for HSB galaxies. At ^B,Um ~ 24, a appears to be 
essentially the same for both LSB and HSB galaxies. A sim- 
ilar pattern can be seen for the L, parameter, except that 
the values for LSB and HSB galaxies appear to be in good 
agreement for fj,B,iim as bright as about 23 mag arcsec"^. 

3.3 Blending 

Since we have used ground-based images of galaxies with ra- 
dial velocities > 840 km s~^ to examine Hll region properties 
in disc galaxies, we cannot ignore the influence blending has 
on the results displayed in Fig. [TT] and 1121 Blending may 
be particularly important when considering the difference 
between the shapes of the Hll region LFs in HSB and LSB 
areas within discs. This is evidenced by the fact that in Fig. 
111! the HSB galaxy LFs have a significantly larger number 
of regions per square kiloparsec than the LSB galaxy LFs 
for fiiim,B < 23, implying that blending may be much more 
important for these LFs. Since the shapes of the LFs for 
HSB and LSB galaxies are roughly similar for /iB,iim^23, 
it is possible that this difference in shape is purely the re- 
sult of the LFs of regions from HSB areas within discs being 
artificially biassing toward higher luminosities by blending. 

The effect of blending can be demonstrated by smooth- 
ing images of more nearby HSB galaxies and reanalysing 
them with Hllphot. For this, we have selected two HSB 
galaxies that were imaged during the October 2000 CTIO 
run, M83 and NGC 1313. These two galaxies were chosen 
because they are relatively nearby {Vr —516 and 457 km 
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Figure 12. For the LFs displayed in Fig. [Sland llll the faint-end 
slope, a (upper), and the value of L* (lower) for the Schechter 
functions fit to the LFs as a function of limiting local B-band 
surface brightness. The parameters from the fits to the LFs using 
all Hll regions are plotted for fiB Um = 0. As in previous figures, 
the results for HSB galaxies are plotted in dark grey and the 
results for LSB galaxies in light grey. 
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Figure 13. For NGC 1313 (upper panels) and M83 (lower panels), the Hll region LFs for four different values of fiB.Hm (see panels) at 
four different simulated radial velocities (see §3.3). The solid lines are Schechter functions fit to the data. 



s~^, respectively) and have ~1,000 Hll regions identified by 
miphot so that LFs can be constructed for each galaxy for 
different values of ^um,B with significant numbers of re- 
gions. In addition to this, M83 and NGC 1313 have two 
distinct morphologies; M83 is a "grand design" spiral galaxy 
(type Sb) and NGC 1313 is an SBd galaxy with a somewhat 
disturbed morphology, especially in Ifa emission. To simu- 
late the effect of blending at different radial velocities, we 
smoothed the B, R, and Ha images of each galaxy with the 
appropriate Gaussian kernel to simulate the size of the see- 
ing disc measured from stars on the narrow-band image at 
10 radial velocites, 560, 600, 700, 800, 900, 1100, 1300, 1500, 
2000, and 2500 km s~^. The images were then re-gridded so 
that the angular size of the pixels relative to that of the see- 
ing disc remained constant. Following this, Gaussian noise 



was added to each image so that the noise in the smoothed 
image was similar to that of the original image. Hllp/iof 
was then used to identify the Hll regions and measure their 
properties on each of these images. 

Following the YlWphot analysis, we constructed Hll re- 
gion LFs for fiB,iim ~0, 22, 23, and 24 for each galaxy at 
each radial velocity. The resulting LFs for Vr =560, 900, 
1500, and 2500 km s"^ are plotted in Fig.[l3]with Schechter 
function fits. In general, as Vr increases, the number of de- 
tected regions decreases and the LFs become more biassed 
toward higher luminosities. To better demonstrate this, we 
have plotted the values of a and L* from the Schechter func- 
tion fits as a function of the physical size of the seeing disc in 
Fig. 1141 For both galaxies, a is not significantly different for 
diflferent values of fiB.Um for ^B,iim < 23. For fiB.Um = 24, 
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a is slightly higher. For M83, a does not appear to change 
significantly with radial velocity while for NGC 1313, a de- 
creases with increasing Vr for jj-Bjim < 23. We therefore 
find no evidence from these two examples that blending has 
caused a to be larger in higher surface brightness areas of 
HSB discs as seen in Fig. [T51 

For both NGC 1313 and M83, L. increases with K for 
HB.iim < 23 and for all values of fiB,iim for M83; L* remains 
roughly constant for NGC 1313 for Ats.iim = 24. For NGC 
1313, the value of L* appears to differ somewhat among 
the different values of i^B,iim with the value for fiB.Um = 24 
being substantially lower at nearly all radial velocities. How- 
ever, the values of L, appear to converge at the lowest radial 
velocities. In contrast, the differences among the values of 
log L* for different limiting surface brightnesses for M83 
are roughly constant with V,-. Linear fits to these data (see 
Fig. I14|) suggest that the values of L, may converge for rel- 
atively small (between ~ 2 and 10 pc) spatial resolutions. 
However, this requires one to extrapolated the trends in the 
lower right panel of Fig. [14] roughly one order of magnitude 
lower than the smallest spatial resolution at Vr — 560 km 
s~^, implying that there is no evidence that blending alone 
may be able to reproduce the trend seen in the lower panel 
of Fig. [T2] for M83. From these two examples, it is clear that 
the values of L, for the higher surface brightness areas of 
HSB galaxies are likely inflated by the effects of blending. 
However, it is also clear that while for some galaxies that 
are similar to NGC 1313, the trend between L, and ^B.Um 
may be entirely caused by blending, for other, more typical 
spiral galaxies like M83, the trend is somewhat unaffected as 
blending appears to increase L, by roughly the same factor 
for all values of ^B,iim- It is therefore likely that the trend 
seen in the lower panel of Fig. [T^Jhas a real physical origin 
and is not entirely the result of blending. 

We note that the Schechter fits were not restricted to 
similar luminosity ranges as those fit to the LFs in Fig. [11] 
because the radial velocity distances for these two galaxies 
are highly uncertain. This has possibly led to the inclusion 
of more faint Hll regions and sharper turn-overs at the faint 
end. We therefore note that the results presented in Fig. 
[14] most likely do not refiect the magnitude of the effect of 
blending on the trend shown in Fig. 1121 However, the results 
are still useful for testing whether or not blending is chiefly 
responsible for the trends seen for HSB galaxies in Fig. 1121 
From the results shown in Fig 1141 it appears that blending 
is likely no the sole cause of these trends. 



4 INTERPRETATION AND DISCUSSION 

As the results presented in Fig. [TT] and [T^] have shown, the 
shape of the Hll region LF appears to be substantially af- 
fected by global processes. In particular, both the faint- 
end slope and L* seem to be roughly independent of lo- 
cal surface brightness for LSB galaxies, while both parame- 
ters seem to decrease as the limiting local B-band surface 
brightness is made fainter for HSB galaxies. In §3.3, we 
demonstrated with two examples that the effects of blend- 
ing have probably inflated the values of L,, possibly to a 
larger degree for higher surface brightness areas. However, 
we also demonstrated that it is likely that blending is not 
the sole cause of the results shown in the lower panel of Fig. 
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Figure 14. For NGC 1313 and M83, the parameters a and L* 
derived from the Schechter functions fit to the LFs measured at 
10 simulated radial velocities (560, 600, 700, 800, 900, 1100, 1300, 
1500, 2000, and 2500 km s'^; see §3.3) for four different values 
of Ms lim versus the physical size of the seeing disc. The dashed 
lines in the bottom panels are weighted, linear least squares fits 
to the data plotted to illustrate the convergence (or lack thereof) 
of the trends between L* and the spatial resolution among the 
different values of fiBdim- 



1121 To analyse the physical origins of the difference in LF 
shapes between LSB and HSB galaxies, we have constructed 
Monte Carlo simula tions of Hll reg i on LF s patterned after 
those developed by lOev fc Clarkel (Il998l) which are simi- 
lar t o those used bylvon Hippelfc Bothun III 19901 ). We note 
that lOev fc Clarke I l| 19981 ) explored different assumptions 
about the time dependence of the Ha luminosity due to 
stars of different masses and about the formation history 
of star clusters. However, they found that assuming a con- 
stant Ha luminosity over the main sequence lifetime of each 
star which abruptly goes to zero and an instantaneous burst 
where all clusters form at once was able to reproduce the 
shapes of typical Hll region LFs. We have therefore adopted 
these two assumptions as well as the same power-laws for 
the distribution of the number of stars per cluster (slope of 
/3 = — 2), the relationships between both Ha luminosity and 
main sequence lifetime and stellar mass (slopes of 5 = 1.5 
and d = —0.7, respective ly), and the initial mass function 
(IMF) of lSalpeterf(|l955h . We have also assumed the same 
upper and lower mass limits for the IMF of 17 and 100 Mq 
and an Ha luminosity and main sequence lifetime for 100 
Mq Stars of 3 x 10^* ergs s"^ and 2.8 Myr. We note that 
these power-law relations break down below the adopted 
lower mass limit for the IMF, which artiflcially introduces a 
turn-over at log L{Ha) ~ 37.7 and we do not consider the 
Monte Carlo LFs below this luminosity. Note that since this 
is a non-physical limit, it was not applied to the fitting of 
Schechter functions to the observed LFs. 

Beyond the assumptions we have made, the two phys- 
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ical factors that determine the shapes of the LPs are the 
age (from birth) of the clusters and the maximum number 
of stars per cluster, N*^max- To explore how these two pa- 
rameters affect the faint-end slope and amount of curvature 
at the bright-end, we have constructed 121 LPs with 1,000 
clusters apiece for 11 ages ranging from to 6 Myr and 11 
values of N^^max ranging from 10 to 1,000. Each simulated 
LP used the same set of random numbers to assign a num- 
ber of stars to each of the 1,000 clusters from the assumed 
power-law distribution, truncated at A^* — Nt,max- We note 
that while ages up to 6 Myr may seem too large for Hll re- 
gions, the instantaneous burst model of lBruzual fc Chariot I 
112003) estimates that at an age of 6 Myr, the flux of ionising 
photons output by a cluster of stars would have dropped by 
a factor of roughly 30, but would be far from negligible. 

An additional factor that must also be considered is 
the uncertainty in the Yia luminosity which varies from 
one galaxy to another and can be rather large because of 
the statistical corrections used for internal dust extinction, 
[Nil] contamination, and stellar absorption [see equation 
(2)]. Also, for galaxies for which the R-band calibration and 
equation (1) were used to calibrate the Hq flux, the uncer- 
tainty in the Ha luminosity is even larger given the uncer- 
tainty in the C parameter in equation (1). This technique 
was used for a larger fraction of the LSB galaxies and may 
therefore contribute to any differences between the shapes of 
the LSB and HSB LPs. To incorporate this into our Monte 
Carlo-simulated LPs, we first computed an estimate of the 
full uncertainty in the Hq calibration for each galaxy. This 
computation included the uncertainties in the parameters 
given in equation (2). It also included the combination of 
the uncertainty in the C parameter in equation (1) and the 
uncertainty in the R-band calibration, or, in the case of the 
majority of the CTIO observations, the uncertainty in the 
calibration solution for the Ha filters derived from obser- 
vations of spectrophotometric standard stars. Pinally, a ra- 
dial velocity uncertainty of 94 km s~^ was assumed for each 
galaxy which is based on the unc ertainty in the correctio n 
for in- fall into the Virgo cluster l|Helmboldt et al. ||2005[ ). 
The resulting uncertainties in log L{Ha) range from about 
0.1 to 0.2 dex. For each group of galaxies, HSB and LSB, the 
fraction of the total number of Hll regions per square kilo- 
parsec among all galaxies in the group was computed for 
each galaxy. This was then used to assign the uncertainty of 
that galaxy to the same fraction of Hll regions in the Monte 
Carlo simulations. This uncertainty was applied to the log 
of the Ha luminosity of each simulated region by adding to 
it a factor of ag where g is a random number drawn from a 
unit normal distribution and a is the assigned uncertainty. 
Since all of the Hll regions within a particular galaxy have 
the same Ha calibration uncertainty, the same value of g was 
used for all simulated regions with the same assigned value 
of a. Each of the 121 simulated LPs was reconstructed twice, 
once using the uncertainties for the LSB group, and once us- 
ing those of the HSB group, resulting in 363 simulated LPs 
in all. 

Schechter functions were fit to all of the simulated LPs, 
including those with and without the modelled effects of 
calibration uncertainty. A sampling of the simulated LPs 
for combinations of six values of star cluster age and six val- 
ues of Nt^max are displayed in Pig. [15] with their Schechter 
fits. Those without the modelled effects of Hq calibration 



uncertainty are plotted in black; those with modelled uncer- 
tainties based on the HSB and LSB groups are plotted in 
grey and light grey, respectively. These plots show that for 
a fixed value of either N^^max or age, the amount of cur- 
vature at the bright-end decreases as the other parameter 
increases. However, it does appear that N^^rnax has a much 
stronger effect on the amount of curvature than the age. 
They also show that the introduction of uncertainties some- 
what reduces the magnitude of this curvature by generally 
flattening the LPs. The magnitude of this effect appears to 
be relatively insignificant for LPs with uncertainties based 
on the data for either HSB or LSB galaxies. 

While other factors such as the assumed IMP, the star 
formation history, and blending also determine the shape 
of the Hll region LP, our Monte Carlo simulations provide 
some insight into the origin(s) of the difference in shapes 
among the LPs plotted in Pig. [11] in terms of age and clus- 
ter size. Prom the simulated LPs shown in Pig. 1151 it appears 
that the power-law shape of the LPs seen throughout LSB 
galaxy discs and in the most LSB regions of HSB discs are 
plausibly the result of clusters forming with relatively large 
values of N^^^max and that the curvature at the bright-end 
seen for the HSB regions of HSB discs requires smaller val- 
ues of Nt,max- However, the paucity of highly luminous Hll 
regions within the LSB regions of dies (see Pig. |9]| suggests 
that it is highly unlikely that N^^max is smaller for clusters 
within HSB regions than for those within LSB regions. Prom 
Pig. 1151 it can be seen that if the clusters within LSB re- 
gions are relatively old (~6 Myr) and the clusters within 
HSB regions are significantly younger (^1 Myr), the value 
of Nf^rnax cau bc roughly the same. Since the effect of blend- 
ing has not been included in the simulated LPs but likely 
affects the observed LPs, we cannot provide accurate esti- 
mates of the typical ages of the clusters within HSB and 
LSB regions using these simulations. However, it appears 
that qualitatively, the simulated LPs suggest that the typi- 
cal star cluster within LSB regions is significantly older than 
the typical HSB region cluster. We note that a similar ex- 
planation has been proposed for the difference in LP shapes 
between arm and inter-arm Hll regions withi n spiral galaxies 
with arm regions being younger on average (|Qev fc Clarke I 
ll998l : lvon Hippel fc Botliun1ll990l ). 

As the LPs in Pig. 1111 demonstrate . the higher surface 
brightness regions of HSB discs tend to have more Hll re- 
gions per square kiloparsec than LSB galaxies have through- 
out their discs. The Monte Carlo-simulated LPs seem to 
imply that these more densely packed areas tend to have 
younger clusters on average. Since the HSB discs have spi- 
ral arms more frequently than the LSB discs in our sam- 
ples, this could point to a fundamental difference between 
spiral-arm driven star formation and the more sporadic star 
formation that is likely prominent in LSB discs, particu- 
larly for Sm and Im galaxies (sec, e.g., Bothun fc Ca l dwelTI 
[Tgsi Ide Blok et aLl 119951 : iHelmboldt et al. ll20oi )7The re- 
sults displayed in Pig. [11] and [12] would further imply that 
there may be a transition from spiral-arm driven star for- 
mation to a more erratic process within HSB discs that typ- 
ically occurs at ~ 24. The galaxy MlOl provides some 
verification of this notion as the luminous Hll regions in the 
outer area of its disc suggest a more erratic process of star 
formation than is s een within the disc w here spiral arms are 
prevalent (see, e.g.. lScowen et aJ~lll992l ). 
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Figure 15. A sampling of the Monte Carlo simulated Hll region LFs for six star cluster ages (marked above the panels) and six values of 
N^^rnax, the maximum number o f stars per cluster (ma rked to the right of the panels). The simulations used 1,000 stellar clusters apiece 
and were pattered after those of lOev fc Clarkel lll998l . see §4.2). The simulated LFs that do not include modelled L{Ha) uncertainties 
are plotted in black. Those that included modelled uncertainties based on the estimated Hq calibration uncertainties (see §4) for LSB 
and HSB galaxies are plotted in light grey and grey, respectively. The plots shown here also include Schcchter function fits [see equation 
(2)]. 



It is also possible that the change in the shape of the 
Hll region LF that appears to occur within HSB galaxies at 
a local surface brightness of /is ~ 24 is linked to a change 
in the cold gas density. For example, the surface brightness 
profile we have measured for M83 reaches = 24 at an ap- 
parent radius of about 4.4 arcmin. From the radial profiles 



presented bv lThilker etHI l|2005l ) , the gas surface density at 
this radius is about 25 Mq pc~^ , which is roughly six times 
larger than the typical critical gas density for sta r formation 
deter mined for similar spiral galaxies (see, e.g., iKennicutt I 
However, at this radius, the profile for the molec- 
ular hydrogen drops quickly, and the slope of log {R) 
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abruptly becomes steeper by a factor of ~ 3. At this radius, 
the gas-phase oxygen abundance of M83 has also dropped by 
factor of ab out two compared t o what it is in the inner parts 
of the disk l|Diaz fc Tosi Ill984l '). This may be the underlying 
cause of the lack of both CO gas and the dust grains needed 
to facilitate the conversion of atomic to molecular hydrogen. 
This example demonstrates that it is plausible that the ap- 
parent transition to a more sporadic mode of star formation 
is linked to a change in physical conditions that has caused 
a sharp decrease in the amount molecular gas that is readily 
available in the outer regions of spiral galaxy discs. 

Overall, the results we have presented demonstrate that 
while local conditions appear to have a large effect on star 
formation within all galaxy discs, it appears that the change 
in Hll region properties from higher to lower density envi- 
ronments is significantly more pronounced within HSB discs 
than it is within LSB discs. The most plausible explana- 
tion appears to be that the star formation within the outer 
parts of HSB discs and throughout LSB discs is more spo- 
radic than that which occurs within the inner regions of 
HSB discs. Since LSB discs tend to have spiral arms less fre- 
quently, little if any molecular gas, and relatively low metal- 
licities, it may be that this sporadic mode of star formation 
tends to dominate within regions that have relatively low 
amounts of dust and molecular gas that are also relatively 
stable against global instabilities like spiral density waves. 
This is also consistent with the physical con ditions of the 
outer regions of spiral ga laxy discs (see, e.g., lOiaz fc Tosi I 
1 1984 iThilker et al~ll2005l ) , which points to a common phys- 
ical origin for the similarly shaped Hll region LFs of both 
LSB discs and the LSB regions of HSB discs. 
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